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Abstract 

These lectures address the dynamics of phase ordering out of equilibrium in 
condensed matter and in quantum field theory in cosmological settings, em- 
phasizing their similarities and differences. In condensed matter we describe 
the phenomenological approach based on the Time Dependent Ginzburg- 
Landau (TDGL) description. After a general discussion of the main exper- 
imental and theoretical features of phase ordering kinetics and the descrip- 
tion of linear (spinodal) instabilities we introduce the scaling hypothesis and 
show how a dynamical correlation length emerges in the large N limit in 
condensed matter systems. The large N approximation is a powerful tool in 
quantum field theory that allows the study of non-perturbative phenomena 
in a consistent manner. We study the exact solution to the dynamics after 
a quench in this limit in Minkowski space time and in radiation dominated 
Friedman-Robertson- Walker Cosmology. There are some remarkable similar- 
ities between these very different settings such as the emergence of a scaling 
regime and of a dynamical correlation length at late times that describe the 
formation and growth of ordered regions. In quantum field theory and cosmol- 
ogy this length scale is constrained by causality and its growth in time is also 
associated with coarsening and the onset of a condensate. We provide a den- 
sity matrix interpretation of the formation of defects and the classicalization 
of quantum fluctuations. 
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I. PHASE ORDERING KINETICS: AN INTERDISCIPLINARY FASCINATING 

PROBLEM 



The dynamics of non-equilibrium phase transitions and the ordering process that occurs 
until the system reaches a broken symmetry equilibrium state play an important role in 
many different areas. Obviously in the physics of binary fluids and ferromagnets (domain 
walls) superffuids (vortex formation), and liquid crystals (many possible defects) to name 
but a few in condensed matter, but also in cosmology and particle physics. In cosmology 
defects produced during Grand Unified Theory (GUT) or the Electro-weak (EW) phase 
transition can act as seeds for the formation of large scale structure and the dynamics of 
phase ordering and formation of ordered regions is at the heart of Kibble's mechanism of 
defect formation [|I]-f|] . Current and future measurements of Cosmic Microwave Background 
anisotropies will determine the nature of the cosmological phase transitions that influenced 
structure formation Also at even lower energies, available with current and forthcoming 
accelerators (RHIC and LHC) the phase transitions predicted by the theory of strong interac- 
tions, Quantum Chromodynamics (QCD) could occur out of equilibrium via the formation 
of coherent condensates of low energy Pions. These conjectured configurations known as 
'Disoriented Chiral Condensates' are similar to the defects expected in liquid crystals or 
ferromagnets and their charge distribution could be an experimental telltale of the QCD 
phase transitions . Whereas the GUT phase transition took place when the Universe was 
about 10 -35 seconds old and the temperature about lO 23 ^, and the EW phase transition 
occured when the Universe was 10~ 12 seconds old and with a temperature 10 15 i^, the QCD 
phase transition took place at about 10~ 5 seconds after the Big Bang, when the temperature 
was a mere 10 12 K. This temperature range will be probed at RHIC and LHC within the 
next very few years. The basic problem of describing the process of phase ordering and the 
competition between different broken symmetry states on the way towards reaching equilib- 
rium is common to all of these situations. The tools, however, are necessarily very different: 
whereas ferromagnets, binary fluids or alloys etc, can be described via a phenomenological 
(stochastic) description, certainly in quantum systems a microscopic formulation must be 
provided. In these lectures we describe a program to include ideas from condensed matter 
to the realm of quantum field theory, to describe phenomena on a range of time and spa- 
tial scales of unprecedented resolution (time scales < 1CT 23 seconds, spatial scales < 1CT 15 
meters) that require a full quantum field theoretical description. 

II. MAIN IDEAS FROM CONDENSED MATTER: 

Before tackling the problem of describing phase ordering kinetics in quantum systems 
starting from a microscopic theory, it proves illuminating to understand a large body of the- 
oretical and experimental work in condensed matter systems [/J- [[TTJ. Although ultimately 
the tools to study the quantum problem will be different, the main physical features to 
describe are basically the same: the formation and evolution of correlated regions separated 
by 'walls' or other structures. Inside these regions an ordered phase exists which eventually 
grows to become macroscopic in size. Before attempting to describe the manner in which 
a given system orders after being cooled through a phase transition an understanding of 
the relevant time scales is required. Two important time scales determine if the transition 
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occurs in or out of equilibrium: the relaxation time of long wavelength fluctuations (since 
these are the ones that order) r re[ (k) and the inverse of the cooling rate t coo \ = T(t)/T(t). 
If r re i(k) << t coo i then these wavelengths are in local thermodynamical equilibrium (LTE), 
but if r re i(k) » t coo i these wavelengths fall out of LTE and freeze out, for these the phase 
transition occurs in a quenched manner. These modes do not have time to adjust locally 
to the temperature change and for them the transition from a high temperature phase to a 
low temperature one occur instantaneously. This description was presented by Zurek [TT| 



analysing the emergence of defect networks after a quenched phase transition. Whereas the 
short wavelength modes are rapidly thermalized (typically by collisions) the long- wavelength 
modes with k « 1/£(T) with £(T) the correlation length (in the disordered phase) become 
critically slowed down. As T — > T c + the long wavelength modes relax very slowly, they fall 
out of LTE and any finite cooling rate causes them to undergo a quenched non-equilibrium 
phase transition. As the system is quenched from T > T c (ordered phase) to T << T c 
(disordered phase) ordering does not occur instantaneously. The length scale of the ordered 
regions grows in time (after some initial transients) as the different broken symmetry phases 
compete to select the final equilibrium state. A dynamical length scale typically emerges 
which is interpreted as the size of the correlated regions, this dynamical correlation length 
grows in time to become macroscopically large [0HTD|- 



The phenomenological description of phase ordering kinetics begins with a coarse grained 
local free energy functional of a (coarse grained) local order parameter M(r) |7],§] which 
determines the equilibrium states. In Ising-like systems this M(r) is the local magnetization 
(averaged over many lattice sites), in binary fluids or alloys it is the local concentration 
difference, in superconductors is the local gap, in superfluids is the condensate fraction etc. 
The typical free energy is (phenomenologically) of the Landau- Ginzburg form: 



F[M}= J d d xl^[VM{x)} 2 + V[M{x)\ 



V[M} = ^r(T)M 2 + ^M 4 ; r(T) = r (T - T c ) (2.1) 

The equilibrium states for T < T c correspond to the broken symmetry states with M = 
±M (T) with 



f for T > T c 

M (T)-j /^ (Tc _ T) i forT<Tc (2-2) 



Below the critical temperature the potential V[M] features a non-convex region with 
d 2 V[M]/dM 2 < for 

- M S (T) < M < M 8 (T) ; M S (T) = ^(T - T c )^ (T < T c ) (2.3) 

this region is called the spinodal region and corresponds to thermodynamically unstable 
states. The lines M S {T) vs. T and M (T) vs. T [see eq. (|2.2j )l are known as the classical 
spinodal and coexistence lines respectively. 

The states between the spinodal and coexistence lines are metastable (in mean-field 
theory). As the system is cooled below T c into the unstable region inside the spinodal, the 
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equilibrium state of the system is a coexistence of phases separated by domains and the 
concentration of phases is determined by the Maxwell construction and the lever rule. 

Question: How to describe the dynamics of the phase transition and the process of 
phase separation? 

Answer: A phenomenological but experimentally succesful description, Time Dependent 
Ginzburg-Landau theory (TDGL) where the basic ingredient is Langevin dynamics @- HIT 



dM{r, t) 
dt 



-T[r,M] 5F \ M \ + v (r,t) 



(2.4) 



with r](f, t) a stochastic noise term, which is typically assumed to be white (uncorrelated) 
and Gaussian and obeying the fluctuation-dissipation theorem: 



(v(r,t)ri(r',t')) = 2TT(r)5 3 (f-r )6(t - t') ; (r ] (f 1 t)} = 



(2.5) 



the averages (■ ■ •) are over the Gaussian distribution function of the noise. There are two 
important cases to distinguish: NCOP: Non-conserved order parameter, with T = To a 
constant independent of space, time and order parameter, and which can be absorbed in a 
rescaling of time. COP: Conserved order parameter with 



-r„ v 



where Tq could depend on the order parameter, but here we will restrict the discussion to 
the case where it is a constant. In this latter case the average over the noise of the Langevin 
equation can be written conservation law 



dM 

~~dT 



j 



-V -J + rj 



d_ 

dt 



( / d 3 fM(r,t)) = 



-r 



6F[M) 
5M 



(2.6) 



where /i is recognized as the chemical potential. Examples of the NCOP are the magnetiza- 
tion in ferromagnets, the gap in superconductors and the condensate density in superfluids 
(the total particle number is conserved but not the condensate fraction), of the COP: the 
concentration difference in binary fluids or alloys. For a quench from T > T c deep into the 
low temperature phase T — > the thermal fluctuations are suppressed after the quench and 
the noise term is irrelevant. In this situation of experimental relevance of a deep quench the 
dynamics is now described by a deterministic equation of motion, 
for NCOP: 



dM 
~dt 



for COP: 



dM 
~dt 



V 2 



5F[M] 
5M 



SF[M] 
' SM 



(2.7) 



(2.8) 



which is known as the Cahn-Hilliard equation In both cases the equations of motion 

are purely diffusive 
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dF_ r 5F[M] dM(f,t) I fd*r(§f) 2 NCOP 

dt ~ J 5M(f,t) dt ~ io | Jd s r (v^) 2 COP { } 

and in both cases ^ < 0. Thus, the energy is always diminishing and there is no possibility 
of increasing the free energy. Thus overbarrier thermal activation cannot be described in 
the absence of thermal noise, which is clear since thermal activation is mediated by large 
thermal fluctuations. The fact that this phenomenological description is purely dissipative 
with an ever diminishing free energy is one of the fundamental differences with the quantum 
field theory description studied in the next sections. 



A. Critical slowing down in NCOP: 



Critical slowing down of long- wavelength fluctuations is built in the TDGL description. 
Consider the case of NCOP and linearize the TDGL equation above the critical temperature 
for small amplitude fluctuations near M — 0. Neglecting the noise term for the moment and 
taking the Fourier transform of the small amplitude fluctuations we find 



dm k (t) 
dt 



k 2 + r (T-T c ) m k (t) 



showing that long-wavelength small amplitud fluctuations relax to equilibrium 
time scale given by 



(2.10) 



on a 



T fc OC 



k 2 + r (T-T c ) 



(2.11) 



As T — > T c + the long-wavelength modes are critically slowed down and relax to equilibrium 
on very long time scales. Therefore a TDGL description leads to the conclusion that if the 
cooling rate is finite, the long-wavelength modes will fall out of LTE and become quenched. 
As the temperature falls below the critical, these modes will become unstable and will grow 
exponentially. 



B. Linear instability analysis: 

Let us consider now the situation for T << T c and neglect the thermal noise. The early 
time evolution after the quench is obtained by linearizing the TDGL equation around a 
homogeneous mean field solution M (t). Writing 

M(f,t) = M (t) + 4= £m fc (t) ^ ( 2 - 12 ) 

where Q is the volume of the system, and considering only the linear term in the fluctuations 
rrikit) the linearized dynamics is the following: COP: for M (t) the conservation gives 

dM {t) 
dt 
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since M Q is the volume integral of the order parameter [see eq.( |2.6| )1 and for the fluctuations 
we obtain 



dm k (t) 
dt 



co(k) rrik(t) ; to(k) = —T k A 



, d 2 V[M) 



k' + 



dM 2 



Mo 



(2.13) 



In the spinodal region 9 v ^ 



dM 2 



Mo 



< there is a band of unstable wave vectors k 2 < 



d 2 V[M] 
dM 2 



Mo 



for which the frequencies are positive and the fluctuations away from the mean 

field grow exponentially. 

NCOP: separate the k ^ from the k = in the linearized equation of motion: 



dM {t) 
dt 

dm k (t) 
dt 



-To 

-r 



dV[M] 



dM 

5F[M] 
5M 



Mo 



m k (t) = -r c 



Mo(t) 



, , d 2 V[M] 



k z + 



dM 2 



Mo 



(2.14) 
(2.15) 



whereas the first equation ( |2.14| ) determines that M Q (t) rolls down the potential hill towards 
the equilibrium solution, the second equation also displays the linear instabilities for the same 
band of wave vectors as in the COP in the spinodal region |M G (t)| < M S (T) [see eq. ( |2.3j )l for 
which the fluctuations grow exponentially in time. Thus in the linearized approximation for 
both NCOP and the COP the spinodal instabilities are manifest as exponentially growing 
fluctuations. These instabilities are the hallmark of the process of phase separation and are 
the early time indications of the formation and growth of correlated regions which will be 
understood in an exactly solvable example below. 



C. The scaling hypothesis: dynamical length scales for ordering 

The process of ordering is described by the system developing ordered regions or domains 
that are separated by walls or other type of defects. The experimental probe to study 
the domain structure and the emergence of long range correlations is the equal time pair 
correlation function 

C(r, t) = (M(f, t)M(0, t)) (2.16) 

where (• • ■) stands for the statistical ensemble average in the initial state (or average over the 
noise in the initial state before the quench) and will become clear(er) below. It is convenient 
to expand the order parameter in Fourier components 

W.*) = 4?£ m *(*) e** 

- 

and to consider the spatial Fourier transform of the pair correlation function 

S(^t) = (m s (t)m_ s (t)> (2.17) 
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known as the structure factor or power spectrum which is experimentally measured by 
neutron (in ferromagnets) or light scattering (in binary fluids) ||12j| . The scaling hypothesis 



introduces a dynamical length scale L(t) that describes the typical scale of a correlated 
region and proposes that 

c{f,t) = f (l^)) ^ S{ ^ t] = Ld{t) 9{kL{t)) (2 ' 18) 

where d is the spatial dimensionality and / and g are scaling functions. Ultimately scaling 
is confirmed by experiments and numerical simulations and theoretically it emerges from a 
renormalization group approach to dynamical critical phenomena which provides a calcula- 
tional framework to extract the scaling functions and the deviations from scaling behavior [|7| . 
This scaling hypothesis describes the process of phase ordering as the formation of ordered 
'domains' or correlated regions of typical spatial size L(t). For NCOP typical growth laws 
are L(t) ~ t 1 / 2 (with some systems showing weak logarithmic corrections) and L(t) ~ t 1 / 3 
for scalar and ~ t 1//4 for vector order parameter in the COP case [J^MIl. 



D. An exactly solvable (and relevant) example: the Large N limit 

We consider the case where the order parameter has ^-components and transforms 
as a vector under rotations in an N-dimensional Euclidean space, i.e. M(f, t) = 
(Mi(r, t), M2(r, i), ■ ■ ■ , M^(r, t)). For N = 1 an example is the Ising model, for N = 2 
superfluids or superconductors (where the components are the real and imaginary part of 
the condensate fraction or the complex gap respectively), N = 3 is the spin one Heisenberg 
antiferromagnet, etc. For N — 1 the topological defects are domain walls (topological in 
one spatial dimension), for N = 2 they are vortices in d = 2 and vortex lines in d = 3, for 
N = d = 3 the topological defects are monopoles or skyrmions which are possible excita- 
tions in Quantum Hall systems and also appear in nematic liquid crystals 0. For iV — > oo 
and fixed d no topological defects exist. However the exact solution of the large TV model 
gives insight and is in fairly good agreement with growth laws for fixed iV systems which 
had been studied experimentally and numerically [7,1C)|. In quantum field theory the non- 



equilibrium dynamics of phase transitions has been studied in Minkowsky and cosmological 



space-times 15-19]. In cosmological space-times it has been implemented to study the col- 



lapse of texture-like configurations [f^ |T3"| , |T4| (see later). The large N limit is an exactly 
solvable limit that serves as a testing ground for establishing the fundamental concepts and 
that can be systematically improved in a consistent 1/N expansion. It provides a consistent 
formulation which is non-perturbative, renormalizable and numerically implementable and 
has recently been invoked in novel studies of non-equilibrium dynamics in quantum spin 
glasses and disordered systems |20| . 



The exact solution for the dynamics in the large N limit, being available both in the 
condensed matter TDGL description of phase ordering kinetics and in Quantum Field The- 
ory in Minkowsky and Cosmological space times, allow us to compare directly the physics 
of phase ordering in these situations. Thus we begin by implementing this scheme in the 
NCOP case for the TDGL description. 
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What is the (• • •) in the equations of the previous section?: consider that before the 
quench the system in in equilibrium in the disordered phase at T >> T c and with a very 
short correlation length (£(T) ~ ^-/T) . The ensemble average in this initial state is therefore 

(M i {f,0)M j {r,0)) = A S ij S 3 (r- r) 

(ilT(f,0)) = (2.19) 

where A specifies the initial correlation. Now consider a critical quench where the system is 
rapidly cooled through the phase transition to almost zero temperature but in the absence 
of explicit symmetry breaking fields (for example a magnetic field). The average of the 
order parameter will remain zero through the process of spinodal decomposition and phase 
ordering. During the initial stages, linear instabilities will grow exponentially with rn l k (i) « 
m l k (0) e^ 1 ; uj(k) = k 2 — r(0) for k 2 < r(0) and at early times 

(mi(t)m j _,(t)} « A e 2uj{k)t (2.20) 

hence fluctuations begin to grow exponentially and eventually will sample the broken sym- 
metry states and the exponential growth must shut-off. The large N limit is implemented 
by writing the potential term in the free energy as 

r(T) A 

V[M] = — \^M 2 + —{M 2 ) 2 - M 2 = M-M (2.21) 

where A is kept finite in the large N limit. We will focus on the NCOP case with a quench 
to zero temperature and rescale the order parameter, time and space as 



M=^i) ; r(0)T t = T ; y^O) x = z (2.22) 
after which the evolution equation for the NCOP case becomes 

(,23) 

where derivatives are now with respect to the rescaled variables. The large N limit is solved 
by implementing a Hartree-like factorization |7[| 

ff — > (if) = N(r] 2 ) no sum over i (2.24) 

Then for each component the NCOP equation becomes 

V 2 + M 2 (t)]r ]i (2.25) 



dr 

M 2 (t) = 1 - ( V 2 ) (2.26) 



the eq. (|2.26|) is a self- consistent condition that must be solved simultaneously with the 
equation of motion for the components. Thus the large iV approximation linearizes the 
problem at the expense of a self-consistent condition. The solution for each component is 
obviously 



S 



lift?) = e- feV+6(T) ; b(r) = f M 2 {r')dr' 

J 



(2.27) 



Consider for a moment that the k = mode is slightly displaced at the initial time, then it 
will roll down the potential hill to a final equilibrium position for which M 2 (oo) 771(00) = (so 
the time derivative vanishes in equilibrium). If 774(00) 7^ is a broken symmetry minimum of 
the free energy, then M 2 (r) — > when r — > 00. This is the statement of Goldstone's theorem 
that guarantees that the perpendicular fluctuations are soft modes. This asymptotic limit 
allows the solution of the self-consistent condition 

M 2 (r) = 1 - (tl?(t)) = 1 - A e 26 ^ / = 1 - A e 26 ^ (8nt)-$ (2.28) 

J (27r) d 

The vanishing of the right hand side in the asymptotic time regime leads to the self-consistent 
solution 



b(r) -^ln 
v ' 4 



M 2 (r) -> A (2.29) 



where r is a constant related to A. This self-consistent solution results in the following 
asymptotic behavior 



%dj)^%m[^J^ T (2-30) 

Introducing the dynamical length scale L(t) = t5 it is straightforward to find the structure 
factor and the pair correlation function 

S(k,t) ocL d (t) e - 2{kL[t))2 (2.31) 



C(r, t) oc e ; L(t) = t* (2.32) 



This behavior should not be interpreted as diffusion, because of the £ d (£) in eqn. ( |2.31| ) 
which is a result of the self-consistent condition. 
Important Features: 

• The 'effective squared mass' M 2 {t) 0: asymptotically there are massless excita- 
tions identified as Goldstone bosons. 

• Since M 2 (t) — > asymptotically, the self-consistent condition results in that (M 2 ) — > 
iW(0)/V% i-e. the fluctuations sample the broken symmetry states, which are equi- 
librium minima of the free energy. These fluctuations begin to grow exponentially at 
early times due to spinodal instabilities. 

• A dynamical correlation length emerges Lit) = t 1 / 2 which determines the size of the 
correlated regions or 'domains'. A scaling solution emerges asymptotically with the 
natural scale determined by the size of the ordered regions. These regions grow with 
this law until they become macroscopically large. Although this a result obtained in 
the large iV limit, similar growth laws had been found for NCOP both analytically 
and numerically for N — 1 etc. M 
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Coarsening: The expression for the structure factor fl2.31|) shows that at large times 
only the very small wavevectors contribute to S(k,t), however the self-consistency 
condition forces the J k d ~ x dk S(k,t) —>■ constant thus asymptotically the structure 
factor is peaked at wavevectors k « L~ l (t) with an amplitude L d (t) thus becoming a 
delta function S(k,t) '^5° 5 d (k). The position of the peak in S(k,t) moving towards 
longer wavelength is the phenomenon of coarsening and is observed via light scattering. 
At long times a zero momentum condensate is formed |T0] and a Bragg peak develops at 
zero momentum, this condensate however grows as a power of time and only becomes 
macroscopic at asymptotically large times. Coarsening is one of the experimental 
hallmarks of the process of phase ordering, revealed for example in light scattering 



12] and is found numerically in many systems 0. Thus the large N limit, although 



not being able to describe topological defects offers a very good description of the 
ordering dynamics. 



III. PHASE ORDERING IN QUANTUM FIELD THEORY I: MINKOWSKI 

SPACE-TIME 

A. A quench in Q.F.T. 

The dynamics is completely determined by the microscopic field theoretical Hamiltonian. 
For a simple scalar theory the Hamiltonian operator is given by 

H = |rf 3 x{^n 2 (f,t) + i[V$(x,t)] 2 + V[$(f,t)]} (3.1) 

where $ is the quantum mechanical field and IT its canonical momentum. We want to 
describe a quenched scenario where the initial state of the system for t < is the ground 
state (or density matrix, see later) of a Hamiltonian for which the potential is convex for 
all values of the field, for example that of an harmonic oscillator, in which case the wave 
function(al) ^[$] is a Gaussian centered at the origin. At t = the potential is changed so 
that for t > it allows for broken symmetry states. This can be achieved for example by 
the following form 

V[<f>] = i m 2 (t)$ 2 + (3.2) 

2m f +m 2 > fort<0 , , 

m(t) = \-m 2 <0 fort>0 (3 ' 3) 

Although in Minkowski space-time this is an ad-hoc choice of a time dependent potential 
that mimics the quench pl| , we will see in the next section that in a cosmological setting 
the mass term naturally depends on time through the temperature dependence and that it 
changes sign below the critical temperature as the Universe cools off. Most of the results 
obtained in Minkowski space-time will translate onto analogous results in a Friedmann- 
Robert son- Walker cosmology. Unlike the phenomenological (but succesful) description of 
the dynamics in condensed matter systems, in a microscopic quantum theory the dynamics 
is completely determined by the Schrodinger equation for the time evolution of the wave 
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function or alternatively the Liouville equation for the evolution of the density matrix in 
the case of mixed states. We will cast our study in terms of a density matrix in general, 
such a density matrix could describe pure or mixed states and obeys the quantum Liouville 
equation 



' at 



H(t),p(t) 



(3.4) 



Question: How does the wave function(al) or the density matrix evolve after a quench? 



B. A simple quantum mechanical picture: 

In order to gain insight into the above question, let us consider a simple case of one 
quantum mechanical degree of freedom q and the quench is described in terms of an harmonic 
oscillator with a time dependent frequency u) 2 {t) = —e(t) u 2 ; oj 2 > with e{t) the sign 
function, so that u 2 (t < 0) > ; u 2 {t > 0) < 0. Furthermore let us focus on the evolution of 
a pure state (the density matrix is simple the product of the wave function and its complex 
conjugate). Consider that at t < the wave function corresponds to the ground state of the 
(upright) harmonic oscillator. For t > the wave function obeys 



1 dt 



Id 2 1 2 2 

2dq~ 2 " 2"° Q 



*[q,t] (3.5) 



Since the initial wave function is a gaussian and under time evolution with a quadratic 
Hamiltonian Gaussians remain Gaussians, the solution of this Schrodinger equation is given 
by 

i&[ q ,t] = N(t) e~^ q2 (3.6) 

« = -La® (3.7) 

Separating the real and imaginary parts of A(t) it is straightforward to find that 
|iV(t)| 4 /ReL4(t)] is constant, a consequence of unitary time evolution. Eq.(|3.8|) can be cast 
in a more familiar form by a simple substitution 

A(t) = =► 4>(t) - u 2 0(t) = (3.8) 



where the equation for was obtained by inserting the above expression for A(t) in ( |3.8|) . 
The solution is <p(t) = a e u ° 1 + b e~^° t featuring exponential growth. This is the quantum 
mechanical analog of the spinodal instabilities described in the previous section. The equal 
time two-point function is given by 

tf)(t)=Al 1 (t) = W)\*™<? uot (3-9) 
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The width of the Gaussian state increases in time (while the amplitude decreases to maintain 
a constant norm) and the quantum fluctuations grow exponentially. As the Gaussian wave 
function spreads out the probability for finding configurations with large amplitude of the 
coordinates increases. These is the quantum mechanical translation of the linear spinodal 
instabilities. When the non-linear contributions to the quantum mechanical potential are 
included the single particle quantum mechanical wave function will simply develop two peaks 
and eventually re-collapse by focusing near the origin undergoing oscillatory motion between 
'collapses' and 'revivals'. In the case of a full quantum field theory there are infinitely many 
degrees of freedom and the energy is transferred between many modes. This simple quantum 
mechanical example paves the way for understanding in a simple manner the main features 
of a quench in the large N limit in quantum field theory, to which we now turn our attention. 

C. Back to the original question: Large N in Q.F.T. 

We now consider the large N limit of a full Q.F.T. in which 

$(£,£) = ($i(z,t),$ 2 (x,t), (3.10) 



and similarly for the canonical momenta IT. The Hamiltonian operator is of the form ( B . 1| ) 
with 

Vm= 1 -m\t)$.$ + ^[$.$] 2 (3.11) 

with m 2 (t) given by ( |3.3|) . Let us focus on the case in which the initial state pure and 
symmetric, i.e. (<&) = 0, with < • • ■ > being the expectation value in this initial state. The 
more complicated case of a mixed state, described by a density matrix is studied in detail 
in [|l6| -18l and the main features are the same as those revealed by the simpler scenario of a 



pure state. The large N limit is implemented in a similar manner as in the TDGL example, 
via a Hartree like factorization 

($• $) 2 -> 2 ($•$)$• I (3.12) 

where the expectation value is in the time evolved quantum state (in the Schrodinger picture) 
or in the initial state of the Heisenberg operators (in the Heisenberg picture). Via this 
factorization the Hamiltonian becomes quadratic at the expense of a self-consistent condition 
as it will be seen below. It is convenient to introduce the spatial Fourier transform of the 
fields as 

£&0 = 4f£*5(*)e** (3-13) 

with Q the spatial volume, and a similar expansion for the canonical momentum U(x,t). 
The Hamiltonian becomes 

h = E {^n r n_ 5 + \w 2 k {t) * s • $_ s } (3.i4) 

k 

W*(t) = m 2 (t) + P + — J — <$ r S_ 5 )(t) (3.15) 
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The problem now has decoupled in a set of infinitely many harmonic oscillators, that are 



only coupled through the self-consistent condition in the frequencies ( 3.15 ). To induce a 
quench, the time dependent mass term has the form proposed in eq. fl3.3p. 

Just as in the simple quantum mechanical case, we consider the initial state to be a 
Gaussian centered at the origin in field space, which is the ground state of the (upright) 
harmonic oscillators for t < 0. Since a Gaussian is always a Gaussian under time evolution 
with a quadratic Hamiltonian, we propose the wave function(al) that describes the (pure) 
quantum mechanical state to be given by 

V[$,t] = U k hf k (t) e-^M-sj ; A k (t = 0) = W k (t < 0) (3.16) 

Time evolution of this wavefunction(al) is determined by the Schrodinger equation: in the 
Schrodinger representation the canonical momentum becomes a differential (functional) op- 
erator, — > —i5/5^_^ and the Schrodinger equation becomes a functional differential 
equation. Comparing the powers of in this differential equation, one obtains the follow- 
ing evolution equations for N k (t) and A k (t) 

j t \nN k {t) = - l -A k (t) (3.17) 

i d -^=Al{t)-Wl{t) (3.18) 

As in the single particle case, the constancy of |iVfc(£)| 4 /Re[v4/ c (t)] is a consequence of unitary 
time evolution. The non-linear equation for the kernel A k (t) can be simplified just as in the 
single particle case by writing 

<Pk{t) 

and taking the expectation value of $ 2 in this state we obtain 

($ r $^)(t)=iV|0 fe (t)| 2 (3.20) 

Hence we find a self-consistent condition much like the one obtained in the large N limit for 
TDGL. The equations for the mode functions and the self-consistent condition for t > are 
therefore given by 

Mt) + [k 2 + M 2 (t)} <f> k {t) = (3.21) 
M\t) = -ml + ±J^\Mt)\ 2 (3.22) 

where the integral in the self-consistent term in ( |3.22|) is simply ($ 2 ). There are two funda- 
mental differences between the quantum dynamics determined by the equations of motion 
and the classical dissipative dynamics of the TDGL phenomeno logical description given in 
sec. II: 



The equations of motion and the self-consistency condition equations ( |3.21| )- (|3T22| ) lead 
immediately to the conservation of energy |15|L6|. 
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• The evolution equations are time reversal invariant. 

These properties must be contrasted to the purely dissipative evolution dictated by the 
TDGL equations as is clear from eq. (2.9). Consider a very weakly coupled theory A << 
1 and very early times, then the self-consistent term can be neglected and we see that 
for k 2 < m\ the modes grow exponentially. This instability again is the manifestation of 
spinodal growth |22|- ^ . p^ ,|l7|] . Since the mode functions grow exponentially, fairly soon, at 
a time scale t s ln(l/A) the self-consistent term begins to cancel the negative mass 

squared and M 2 (t) becomes smaller. We find numerically that this effective mass vanishes 
asymptotically, as shown in Fig. 1. 



D. Emergence of condensates and classicality: 

The physical mechanism here is similar to that in the classical TDGL, but in terms of 
quantum fluctuations. The quantum fluctuations with wave vectors inside the spinodally 
unstable band grow exponentially, these make the ($ 2 ) self-consistent field to grow non- 
perturbatively large until when ($ 2 ) ps rn^/X when the self-consistent (mean) field begins to 
be of the same order as m\ (the tree level mass term). At this point the quantum fluctuations 
become non-perturbatively large and sample field configurations near the equilibrium min- 
ima of the potential. The spinodal instabilities are shutting off since the effective squared 
mass M 2 (t) is vanishing. 

When M 2 (t) vanishes, the equations for the mode functions become those of a free 
massless field, with solutions of the form (j>k{t) = Aj, e lkt + e~ tkt , whereas for the k = 
mode the solution must be of the form <fio{t) = a + bt with a; b ^ since the Wronskian of 
the mode function and its complex conjugate is a constant. This in turn determines that 
the low k (long wavelength) behavior of the mode functions is given by 

sin kt 

(pk\t) = a cos kt + b — - — (3.23) 

K 

This behavior at long wavelength has a remarkable consequence: at very long time the 
power spectrum \<pk{t)\ 2 , which is the equivalent of S(k,t) for TDGL (see eq. ( |2.17| )) is 



dominated by the small fc-region, in particular k « 1/t, with an amplitude that grows 
quadratically with time. Then the structure factor S(k,t) = \4>k(t)\ 2 features a peak that 
moves towards longer wavelengths at longer times and whose amplitude grows with time in 
such a way that asymptotically J °° k 2 S(k,t)dk/2n 2 — > tUq/X and the integral is dominated 
by a very small region in k that gets narrower at longer times. This is the equivalent of 
coarsening in the TDGL solution in the large N limit, where the asymptotic time regime 
was dominated by the formation of a long-wavelength condensate. Fig. 2 shows the power 
spectrum at two (large) times displaying clearly the phenomenon of coarsening and the 
formation of a non-perturbative condensate. 

The pair correlation function can now be calculated using this power spectrum JTT 



C(r,t) = ^r k sin kr |0 2 (t)| dk . (3.24) 
2 7i z r Jo 
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At long times and distances the integral is dominated by the very long wavelength modes, 
in particular by the term oc sin[kt]/k of 4>k{t), hence the integral can be done analytically 
and we find 

C(f,t) = -e(2t-r) (3.25) 
r 

with A a constant. This is a remarkable result: the correlation falls off as 1/r inside domains 
that grow at the speed of light. This correlation function is shown in Fig. 3 at several 
different (large) times. This correlation function is of the scaling form: introducing the 
dynamical length scale L(t) = t it is clear that |l7j 



C{r,t)KL-\t)f{r/L{t)) ; f(s) = 9(2 s) (3.26) 

s 

We interpret these 'domains' as being a non-perturbative condensate of Goldstone 
bosons, with a non-perturbatively large number of them oc 1/A, such that the mean square 
root fluctuation of the field samples the (non-perturbative) equilibrium minima of the poten- 
tial. In particular an important conclusion of this analysis is that the long-wavelength modes 
acquire very large amplitudes, their phases vary slowly as a function of time (for k << 1/t), 
therefore these fluctuations which began their evolution as being quantum mechanical, now 
have become classical. 



E. O.K. ..O.K. but where are the defects? 

At this point our analysis begs this question. To understand the answer it is convenient 
to back track the analysis to the beginning. The initial quantum state is given by a the 
wave-function(al) (|3.16|) , thus the most probable field configurations found in this ensemble 
are those whose spatial Fourier transform are given by 

\$ k \ oc ^=^= oc ^=^= (3.27) 



W k (t<0) Jk 2 + 



m. 



(restoring % would multiply by Vh). Then typical long- wavelength field configurations 
that are represented in the quantum ensemble described by this initial wave-function(al) are 
of rather small amplitude. The initial correlations are also rather short ranged on scales 
tuq 1 . Under time evolution the probability distribution is given by 



v[$,t] = |^[<M| 2 = n^n, ||iv fc (t)| 2 e (3.28) 

At times longer than the regime dominated by the exponential growth of the spinodally 
unstable modes, the power spectrum |0 2 (t)| 2 obtains the largest support for long wavelengths 
k « ttLq and with amplitudes ~ m^/A. Therefore field configurations with typical spatial 
Fourier transform 4>k{t) are very likely to be found in the ensemble. These field configurations 
are primarily made of long-wavelength modes and their amplitudes are non-perturbatively 
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large, of the order of the amplitude of the fields in the broken symmetry minima. A typical 
such configuration can be written as 



&(Z,t)t yP icai « E \Mt)\ cos[k -x + 5i] (3.29) 

k 

where the phases 5| are randomly distributed with a Gaussian probability distribution since 
the density matrix is gaussian in this approximation. We note that a particular choice of 
these phases leads to a realization of a likely configuration in the ensemble that breaks trans- 
lational invariance. In fact translations can be absorbed by a change in the phases, thus 
averaging over these random phases restores translational invariance. Since the quantum 
state (or density matrix) is translational invariant a particular spatial profile for a field con- 
figuration corresponds to a particular representative of the ensemble. Combining all of the 
above results together we can present the following consistent interpretation of the ordering 
process and the formation of coherent non-perturbative structures during the dynamics of 
symmetry breaking in the large N limit [jnj : 

• The early time evolution occurs via the exponential growth of spinodally unstable 
long wavelength modes. This unstable growth leads to a rapid growth of fluctuations 
($ 2 )(t) which in turn increases the self-consistent contribution and tends to cancel the 
negative mass squared. The effective mass of the excitations —ttIq + -^(^ 2 )(t) — ► 
and the asymptotic excitations are Goldstone bosons. 

• At times larger than the spinodal time t s ~ uiq 1 ln(l/X), the effective mass vanishes 
and the power spectrum or structure factor S(k,t) = \4>k{t)\ 2 displays the features of 
coarsening: a peak that moves towards longer wavelengths and increases in amplitude, 
resulting in a long-wavelength condensate at asymptotically long times. 

• For large time a dynamical correlation length emerges L(t) = t and at long distances 
the pair correlation function is of the scaling form C(r,t) oc L^ 1 (t) f (r / L(t)) . The 
length scale L(t) determines the size of the correlated regions and determines that 
these regions grow at the speed of light. Inside these regions there is a non-perturbative 
condensate of Goldstone bosons with a typical amplitude of the order of the value of 
the homogeneous field at the equilibrium broken symmetry minima. 

The similarity between these results and those of the more phenomenological TDGL descrip- 
tion in condensed matter systems is rather striking. The features that are determined by 



the structure of the quantum field theory are ||17|| : i) the scaling variable s = r/t with equal 



powers of distance and time is a consequence of the Lorentz invariance of the underlying 
theory, ii) the fact that the pair correlation function vanishes for r > It is manifestly a con- 
sequence of causality. An analysis of the correlations and defect density during the spinodal 
time scale has been performed in [p5l and related recent studies had been performed in |26| . 



IV. PHASE ORDERING IN QUANTUM FIELD THEORY II: FRW COSMOLOGY 
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A. Cosmology 101 (the basics): 



On large scales > 100 Mpc the Universe appears to be homogeneous and isotropic as 
revealed by the isotropy and homogeneity of the cosmic microwave background and some of 
the recent large scale surveys ||. The cosmological principle leads to a simple form of the 
metric of space time, the Friedmann-Robertson- Walker (FRW) metric in terms of a scale 
factor that determines the Hubble flow and the curvature of spatial sections. Observations 
seem to favor a flat Universe for which the space time metric is rather simple: 

ds 2 = dt 2 -a 2 {t)dx 2 (4.1) 

the time and spatial variables t, x in the above metric are called comoving time and spatial 
distance respectively and have the interpretation of being the time and distance measured 
by an observer locally at rest with respect to the Hubble flow. At this point we must note 
that physical distances are given by l p hy S (t) = a if) An important concept is that of 
causal (particle) horizons: events that cannot be connected by a light signal are causally 
disconnected. Since light travels on null geodesies ds 2 = the maximum physical distance 
that can be reached by a light signal at time t is given by 

r* dt' 

d H (t) = a(t) / — - (4.2) 
Jo a(t') 

It will prove convenient to change coordinates to conformal time by defining a conformal 
time variable 

r* dt' 

V= ~-^ds 2 = C 2 ( V ) {drf - dx 2 ) ■ C( V ) = a(t(rj)) (4.3) 
Jo a{v) 

in terms of which the causal horizon is simply given by dn{v) — C(rj) rj and physical 
distances as x p h ys = C(n) x. This metric is of the same form as that of Minkowski space 
time. For energies well below the Planck scale Mp\ « 10 19 Gev gravitation is well described 
by classical General Relativity and the Einstein equations: 

BT - - g^R = -4- (4.4) 

where we have been cavalier and set c = 1 (as well as h — 1). R^ u is the Ricci tensor, R 
the Ricci scalar and T^ u the matter field energy momentum tensor. The above equation 
is classical but one seeks to understand the dynamics of the Early Universe in terms of a 
quantum field theory that describes particle physics, thus the question: what is exactly the 
energy momentum tensor?, in Einstein's equations it is a classical object, but in QFT it is 
an operator. The answer to this question is: gravity is classical, fields are quantum mechan- 
ical, but T^ u — > (T Ml/ ), i.e. it is the expectation value of a quantum mechanical operator 
in a quantum mechanical state. This quantum mechanical state, either pure or mixed is 
described by a wave-function(al) or a density matrix whose time evolution is dictated by the 
quantum equations of motion: the Schrodinger equation for the wave functions or the quan- 
tum Liouville equation for a density matrix. Consistency with the postulate of homogeneity 
and isotropy requires that the expectation value of the energy momentum tensor must have 
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the fluid form and in the rest frame of the fluid takes the form (T 7 ^) = diagonal(p,p, p,p) 
with p the energy density and p the pressure. The time and spatial components of Einstein's 
equations lead to the Friedman equation 

a 2 (t) 8n /A 

W) = ™i/' (t) (4 ' 5) 

2^ + ^l = -^PW (4-6) 



a(t) a?(t) M\ 



2 

PI 



Combining these two equations one arrives at a simple and intuitive equation which is 
reminiscent of the first law of thermodynamics: 

^(pa 3 (t)) = -p^- =► p + 3% + p) = (4.7) 
at at a 

The alternative form shown on the right hand side of ( |4.7|) is the covariant conservation 
of energy. Since the physical volume of space is V a 3 (t) (with V the co moving volume) 
the above equation is recognized as dU = —p dV which is the first law of thermodynamics 
for adiabatic processes. To close the set of equations and obtain the dynamics we need 
an equation of state p = p(p): two very relevant cases are: i) radiation dominated (RD) 
with p = p/3 and matter dominated (MD) p = (dust) Universes. In our study we will 
focus on the RD case. The equation of state for RD is that for blackbody radiation for 
which the entropy is 5* = CVT 3 (with C a constant). Since V(t) = Vq a 3 (t) is the physical 
volume, the equation (^4.7|) which dictates adiabatic (isoentropic) expansion leads to a time 
dependence of the temperature: T(t) = T /a(t). Now the cooling is done by the expansion 
of the Universe and a phase transition will occur when the Universe cools below the critical 
temperature for a given theory. For the GUT transition T c « 10 16 Gev « 10 29 i^, for the 
EW transition T c w 100 Gev w 10 15 K. Returning now back to the large iV study of the 
dynamics of phase transitions, we can include the effect of cooling by the expansion of the 
Universe by replacing the time dependent mass term m 2 (t) in fl3.11p by 



m 2 (t) = ml 



T\t) 

c 



T « = J) < 4 ' 8 ' 



This form is consistent with the Landau- Ginzburg description including the time depen- 
dence of the temperature via the isentropic expansion of the Universe, but perhaps more 
importantly it can be proven in a detailed manner from the self-consistent renormalization 
of the mass in an expanding Universe |15| . Thus the large N limit in a RD FRW cosmology 
will be studied by using the potential ( [3.1 1| ) but with the time dependent mass given by 



B. Large N in a RD FRW Cosmology 

The large N limit is again implemented via the Hartree-like factorization ( |3.12| ) perform- 
ing the spatial Fourier transforms of the fields and their canonical momenta and including 
the proper scale factors, the Hamiltonian now becomes 
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H{t) = E {^y % • n_ fe - + w 2 k {t) $ z • $_ fc - 



k 

k 2 



A 



(4.9) 
(4.10) 



where now the expectation value is in terms of a density matrix p[$(7), $(T); i] since we are 
considering the case of a thermal ensemble as the initial state. 

We propose the following Gaussian ansatz for the functional density matrix elements in 



the Schrodinger representation [15 



A k (t) 



(4.11) 



This form of the density matrix is dictated by the hermiticity condition pt[$, $, £] = 
p*[$,$,£]; as a result of this, Bkit) is real. The kernel Bkit) determines the amount of 
mixing in the density matrix, since if B k = 0, the density matrix corresponds to a pure state 
because it is a wave functional times its complex conjugate. The kernels Aj.(0) ; Bj,(0) are 
chosen such that the initial density matrix is thermal with a temperature Tj > T c |15[ . Fol- 
lowing the same steps as in Minkowski space time, the time evolution of this density matrix 
can be found in terms of a set of mode functions <fik{t) that obey the following equations of 
motion and self-consistency condition 



<f> k {t) + 3 - <p k (t) + 
a 



M 2 {t) 



T? 



T 2 a 2 (t) 



a?(t) 
- 1 



+ M 2 (t) 
A 

+ 2 



<Pk{t) = 



(2tt 



(4.12) 
(4.13) 



This equations can be cast in a more familiar form by changing coordinates to conformal 
time (see eq. (|4.3|) ) and (conformally) rescaling the mode functions <j>k(t) = fk{xi)/C{rj) to 
obtain the following equations for the conformal time mode functions fk{v) in an RD FRW 
cosmology 



fk(v)+ k 2 + C 2 ( V )M 2 ( V ) f k ( V )=0 



(4.14) 



where primes now refer to derivatives with respect to conformal time. For RD FRW C(rj) = 
l + ?7/2 (in units of ttiq 1 which is the only dimensionful variable). The above equations of 
motion has now an analogous form as those solved in the case of Minkowski space-time. 

As the temperature falls below the critical the effective squared mass term becomes 
negative and spinodal instabilities trigger the process of phase ordering. This results in that 
the quantum fluctuations quantified by (<5 2 ) grow exponentially. These spinodal instabilities 
make the self-consistent field grow at early times and tends to overcome the negative sign 
of the squared mass, eventually reaching an asymptotic regime in which the total effective 
mass M 2 (rj) vanishes. 

Again this behavior determines that the fluctuations are sampling the equilibrium broken 
symmetry minima of the initial potential, i.e. ($ 2 ) 



2Nml 
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Although, just as in Minkowski space-time the effective mass vanishes asymptotically, the 



non-equilibrium evolution is rather different. We find numerically [19fl that asymptotically 
the effective mass term C 2 (r])M 2 (r]) vanishes as — 15/Ar] 2 . 

Fig. 4 displays C 2 (f])M 2 (i]) as a function of conformal time for the case of Tj/T c = 1.1 
with T c ocm /v4 

We see that at very early time the mass is positive, reflecting the fact that the initial 
state is in equilibrium at an initial temperature larger than the critical. As time evolves the 
temperature is red-shifted and cools and at some point the phase transition occurs, when 
the mass vanishes and becomes negative. 

Figure 5 displays 21 ^ m i (<& 2 )(r]) vs. r\ in units of wig 1 for |f = 3, g = 10~ 5 for an R.D. 
Universe. Clearly at large times the non-equilibrium fluctuations probe the broken symmetry 
states. 

This particular asymptotic behavior of the mass determines that the mode functions 
fk(v) grow as rfl 2 for k < 1/rj and oscillate in the form e ±lkv for k > l/rj. This behavior is 
confirmed numerically [|I9[| . We find both analytically and numerically that asymptotically 
the mode functions are of the scaling form 

AM^W (415) 

Where A is a numerical constant and ^(a?) is a Bessel function. 

Figure 6 displays i]~ 5 \fk{v)\ 2 as a function of the scaling variable kr\ revealing the scaling 
behavior. 

It is remarkable that this is exactly the same scaling solution found in the classical non- 
linear sigma model in the large N limit and that describes the collapse of textures [O], and 



also within the context of TDGL equations in the large N limit applied to cosmology |T4 . 

The growth of the long-wavelength modes and the oscillatory behavior of the short wave- 
length modes again results in that the peak of the structure factor S(k,r]) = \fk(v)\ 2 moves 
towards longer wavelengths and the maximum amplitude increases. This is the equivalent 
of coarsening and the onset of a condensate. 

Although quantitatively different from Minkowsky space time, the qualitative features 
are similar. Asymptotically the non-equilibrium dynamics results in the formation of a non- 
perturbative condensate of long-wavelength Goldstone bosons. We can now compute the 
pair correlation function C (r, 77) from the mode functions (|4.15|) and find that it is cutoff 
by causality at r = 2rj. The correlation function is depicted in Fig. 4 for two different 
(conformal) times. 

The scaling form of the pair correlation function is 

C(r, 77) = rf x(r/v) 

where x( x ) * s a hump-shaped function as shown in fig. 7. 

Clearly a dynamical length scale L{rj) = i] emerges as a consequence of causality, much 
in the same manner as in Minkowsky space time. The physical dynamical correlation length 
is therefore given by Cphysiv) = C(r])L(r]) = d H (t), that is the correlated domains grow again 
at the speed of light and their size is given by the causal horizon. The interpretation of this 
phenomenon is that within one causal horizon there is one correlated domain, inside which 
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the mean square root fluctuation of the field is approximately the value of the equilibrium 
minima of the tree level potential, this is clearly consistent with Kibble's original observation 
P]H]. Inside this domain there is a non-perturbative condensate of Goldstone bosons ||19|| . 



V. CONCLUSIONS AND LOOKING AHEAD 

In these lectures we have discussed the multidisciplinary nature of the problem of phase 
ordering kinetics and non-equilibrium aspects of symmetry breaking. Main ideas from con- 
densed matter were discussed and presented in a simple but hopefully illuminating frame- 
work and applied to the rather different realm of phase transitions in quantum field theory as 
needed to understand cosmology and particle physics. The large iV approximation has pro- 
vided a bridge that allows to cross from one field to another and borrow many of the ideas 
that had been tested both theoretically and experimentally in condensed matter physics. 
There are, however, major differences between the condensed matter and particle physics- 
cosmology applications that require a very careful treatment of the quantum field theory 
that cannot be replaced by simple arguments. The large N approximation in field theory 
provides a robust, consistent non-perturbative framework that allows the study of phase 
ordering kinetics and dynamics of symmetry breaking in a controlled and consistently im- 
plementable framework, it is renormalizable, respects all symmetries and can be improved 
in a well defined manner. This scheme extracts cleanly the non-perturbative behavior, the 
quantum to classical transition and allows to quantify in a well defined manner the emer- 
gence of classical stochastic behavior arising from non-perturbative physics. The emergence 
of scaling and a dynamical correlation length are robust features of the dynamics and the 
Kibble-Zurek scenario describes fairly well the general features of the dynamics, albeit the 
details require careful study, both analytically and numerically. 

Of course this is just the beginning, we expect a wealth of important phenomena to be 
revealed beyond the large N, such as the approach to equilibrium, the emergence of other 
time scales associated with a hydrodynamic description of the evolution at late times and 
a more careful understanding of the reheating process and its influence on cosmological 
observables. Although within very few years the wealth of observational data will provide 
a more clear picture of the cosmological fluctuations, it is clear that the program that 
pursues a fundamental understanding of the underlying physical mechanisms will continue 
seeking to provide a consistent microscopic description of the dynamics of cosmological phase 
transitions. 
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FIG. 1. M 2 {t) vs. t, g = 1(T 
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FIG. 3. grC{r,r) vs r/\m \ for t/\m \ = 100, 200, 300 for g = 10~ 7 . 
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FIG. 5. 2A^(^ 2 )( r ?) vs - V (conformal time in units of ttiq 1 ) for ^ = 3, g = 10 5 . R.D. 
Universe. 
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FIG. 7. C(r,7]) vs. r for rj = 250,400 (in units of m^ 1 ) for J = 3, g = 10~ 5 . R.D. FEW 
Universe. 
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